Singlet-triplet dephasing in asymmetric quantum dot molecules 
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We discuss pure dephasing of singlet-triplet superpositions in two-electron double quantum dots 
due to elastic phonon scattering. We generalize our previous results to a system built of two non- 
identical dots. We show that the asymmetry must be very strong in order to considerably affect the 
dephasing rate. 
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Dephasing of singlet-triplet superpositions in two- 
electron quantum dot molecules is important for a pos- 
sible implementation of quantum information processing 
in semiconductor systems [lj. In a recent work [2|, we 
showed that elastic phonon scattering via virtual tran- 
sitions to doubly occupied states, which is only possible 
in a singlet configuration, induces distinguishability of 
spin configurations and, therefore, leads to pure dephas- 
ing of spin superpositions. In Ref. 0, we studied a sys- 
tem of two identical dots. However, as quantum dots are 
artificial systems, one has to take into account unavoid- 
able inhomogeneity of dot parameters when modeling the 
properties of the system. Therefore, in the present con- 
tribution, we generalize our previous result and study 
the phonon-induced dephasing process in an asymmetric 
QDM. We show that in the asymmetric QDM, an ad- 
ditional dephasing channel appears, as compared to the 
symmetric one. Nonetheless, the dephasing rate is very 
weakly affected by the asymmetry unless the latter be- 
comes very strong. 

The system under consideration is composed of two 
electrons in an asymmetric quantum dot molecule 
(QDM) built from two different gate-defined quantum 
dots @, |3] . The electrons are coupled to phonons by 
the usual charge-phonon interactions (deformation po- 
tential and piezoelectric couplings). We do not take any 
spin-environment interactions into account (neither di- 
rect, with nuclear spins, nor indirect, via spin-orbit cou- 
pling) . The Hamiltonian of the system is then 

H = HbQD + H p h + -flint. 

The first term describes the electrons and has the form 
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respectively, and s =|,J. labeling the spin orientation. 
The first term in Eq. {I} accounts for the energy differ- 
ence between single-electron states in the two dots. The 
second term represents single-particle inter-dot tunnel- 
ing. The third term describes the Coulomb interaction, 
with Vijki = Vjuk = Vkuj = Vikji (the wave functions 
may be chosen such that the matrix elements are real). 

The Hamiltonian of the phonon reservoir is given by 
#ph = J2k,x ftu k ,\b k > bk,\, where b h ,\,b\ x are annihila- 
tion and creation operators for a phonon from a branch 
A with a wave vector k and Tiujk,\ are the corresponding 
energies. The electron-phonon interaction is described 
by 



S,2 fc,A 



where F^ R (k) are coupling constants. We use the usual 
coupling constants for confined charges (see Refs. 0, HI 
for explicit expressions). 

The Hamiltonian i?DQD can be diagonalized in a sim- 
ple way. The resulting singlet eigenstates can be written 
in the form 
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where ai S , a\ s are the electron annihilation and creation 
operators with i = L, R denoting the left and right dot, 



where |(±)5> = [|(2,0)5) ± |(0, 2)S)]/y/2 and \{m,n)S) 
denotes the configuration with m electrons in the left 
dot and n electrons in the right one. The corresponding 
energies are denoted by Eq, Ei,E + . 

Here 9 describes mixing of the |(1, 1)5) and |(+)5) con- 
figurations due to tunneling and Coulomb interactions. 
In the present work we fix 6 = 0.5, which corresponds to 
t 1 /(E 1 - E ) = 0.2. 

The angle <f> accounts for state mixing due to system 
asymmetry. The latter may be due to a difference be- 
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tween the energies of doubly occupied configurations in 
the two dots or to asymmetry of off-diagonal Coulomb 
elements. Accordingly, we define two parameters 
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Neglecting the exchange terms compared to direct 
Coulomb energies, one finds sin0 w [sin(0/2)A — 
cos(6/2)r ] ]/(E 1 -E ). 

The carrier-phonon coupling for the two-electron sys- 
tem is written in the eigenbasis of -Hdqd , 

H int = FS(k)\a)(p\(b k , x + bl k 



-fc.A ) ' 



with a, (3 — So, S\, S+. Then, the phonon spectral den- 
sities involving the low-energy state \Sq) are found, 

fc,A 

where a,[3 = S\,S+ and n(u>) is the Bose distribution. 
In contrast to thy previously studied symmetric case, two 
channels of phonon scattering are now present, corre- 
sponding to the two different excited states (Si and S+) 
through which the phonons can scatter. The interference 
of different scattering paths is reflected by the presence of 
the whole family of spectral densities out of which only 
one, Rs 1 Si( UJ ) survives for <j> — > and reduces to the 
spectral density found in Ref. Q . 

Following the method worked out in our previous pa- 
per 0] , we calculate the dephasing rate using the 4th or- 
der time-convolutionless equation for the evolution of the 
density matrix describing the state of the two-electron 
system. From this, we extract the Markov limit, which 
yields the dephasing rate due to two-phonon (scatter- 
ing) processes. Similarly as in the symmetric case 
this process dominates over the one-phonon real transi- 
tions at low temperatures. The resulting dephasing rate 
is given by 
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where uj a = (E a — Eq)/% and V denotes the Cauchy 
principal value. The above equation generalizes the pre- 
viously found formula to the asymmetric case. 

In the calculations, we set the energy differences Ei — 
E Q = 0.9 meV, E + - E Q = 1.0 meV, the inter-dot dis- 
tance D = 300 nm, and use the material parameters as 
in Ref. 0. 



The pure dephasing rates resulting from the scattering 
process are shown in Figs. [2a, c) as a function of A for a 
1 
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FIG. 1: Dephasing rate: (a,c) as a function of A for 77 = 
(solid line), 0.2 meV (dashed line), 0.3 meV (dotted line), and 
0.4 meV (dash-dotted line); (b,d) as a function of 77 for A = 
(solid line), 0.2 meV (dashed line), 0.4 meV (dotted line), and 
0.6 meV (dash-dotted line). In (a) and (b), T = 0.5 K; in (c) 
and (d), T — IK. 



few values of 77 and in Fig. []](b,d) as a function of 77 for 
a few values of A. In all the cases, the dephasing time 
varies from about 1 /is to several /zs. As long as the asym- 
metry is small, the dephasing rate is rather insensitive to 
it. Only for strong asymmetry, the dephasing rate grows 
considerably One should keep in mind that the largest 
values of 77 shown in Fig. [1] are rather unrealistic since 
this parameter describes the variation of exchange-like 
Coulomb terms which are themselves small. In a real 
structure, 77 is likely to remain within a tenth of meV. 
It is clear from Fig. [T] that in this parameter range the 
dephasing rate is almost completely insensitive to 77, as 
well as to A. 
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Foundation (Grant no. 202/07/J051). 



[1] 

[2] 
[3] 



Electronic address: pawel.machnikowski@pwr.wroc.pl 
D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 

(1998). 

K. Roszak and P. Machnikowski. larXiv:0903. 07831 

J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. 

Yacoby, M. D. Lukin, CM. Marcus, M. P. Hanson, and 

A. C. Gossard, Science 309, 2180 (2005). 

F. H. L. Koppens, J. A. Folk, J. M. Elzerman, R. Hanson, 

L. H. Willems van Beveren, I. T. Vink, H. P. Tranitz, W. 

Wegscheider, L. P. Kouwenhoven, and L. M. K. Vander- 

sypen, Science 309, 1346 (2005). 

A. Grodecka, P. Machnikowski, and J. Forstner, Phys. 
Rev. B 78, 085302 (2008). 



